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ABSTRACT 

Perturbations of rotating relativistic stars can be classified by their behavior 
under parity. For axial perturbations (r-modes), initial data with negative 
canonical energy is found with angular dependence e* m ^ for all values of m > 2 
and for arbitrarily slow rotation. This implies instability (or marginal stability) 
of such perturbations for rotating perfect fluids. This low m-instability is 
strikingly different from the instability to polar perturbations, which sets in first 
for large values of m. The timescale for the axial instability appears, for small 
angular velocity fl, to be proportional to a high power of fh As in the case of 
polar modes, viscosity will again presumably enforce stability except for hot, 
rapidly rotating neutron stars. 

This work complements Andersson’s numerical investigation of axial modes 
in slowly rotating stars. 


Subject headings: instabilities — relativity — stars: oscillations — stars: 
rotation 


1. Introduction 

A rotating star is invariant under reflection in the equatorial plane and rotations about 
the axis of symmetry. It is therefore invariant under parity, the composition of reflection 
in the equatorial plane and rotation by 7r, and perturbations of a rotating star can be 
chosen to have definite parity. That is, the space of solutions to the Einstein-perfect-fluid 
equations, linearized about the equilibrium star, is a direct sum of dynamically preserved 
subspaces with parity eigenvalues ±1. [] 


1 It would be natural to call these “odd” or “even” perturbations, but some of the literature in 
relativistic astrophysics adopts the Regge-Wheeler (1957) terminology, using the term “odd-parity” not 
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Perturbations of a spherical star can be divided into two classes, axial and polar, 
depending on their behavior under parity. Where polar tensor fields on a 2-sphere can 
be constructed from the scalars Y) m , their gradients V aY f n , and the metric cab cm the 
2-sphere, axial fields involve the pseudo-tensor €ab, and their behavior under parity is 
opposite to that of Y t m . Thus axial perturbations of odd l are invariant under parity, and 
axial perturbations with even l change sign. If a mode varies continuously along a sequence 
of equilibrium configurations that starts with a spherical star and continues along a path 
of increasing rotation, the mode will be called axial if it is axial for the spherical star. 

Its parity cannot change along the sequence, but l is well-defined only for modes of the 
spherical configuration. 

Axial perturbations of a spherical stellar model are time-independent convective 
currents that do not change the density and pressure of the star and do not couple to 
gravitational waves. Along a continuous sequence of relativistic perfect-fluid models, a mode 
changes from stable to unstable when its frequency is zero; because they have vanishing 
frequency in nonrotating models, axial modes could thus be unstable for any model with 
nonzero rotation. Andersson’s numerical work and the computation we present here shows 
that this is the case: Axial modes for every value of m are unstable for arbitrarily slow 
rotation. 

In spherical stars, gravitational radiation removes positive angular momentum L z 
from a mode moving in the positive (j) direction and negative angular momentum from 
a backward-moving mode; and it therefore damps all time-dependent nonaxisymmetric 
modes. Once the angular velocity of the star is sufficiently large, however, a mode that 
moves backward relative to the star is dragged forward relative to an inertial observer. 
Gravitational radiation will then remove positive angular momentum from the mode. 

But a mode that moves backwards relative to the fluid has negative angular momentum, 
because the perturbed fluid does not rotate as fast it did without the perturbation. The 
radiation thus removes positive angular momentum from a mode whose angular momentum 
is negative. By making the angular momentum of the perturbation increasingly negative, 
gravitational radiation drives the mode (Friedman and Schutz 1978, Friedman 1978, 
Chandrasekhar 1970). 

For slowly rotating Newtonian stars, the frequencies of axial modes are known (see 
Sect. |2|), and they lie in this unstable regime. In general relativity, Andersson (1997) 
numerically computes the real part of the frequency of these modes in a slow-rotation 


for a perturbation that changes sign under parity but for an axial perturbation, whose behavior is opposite 
to that of Y t m . 
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approximation, while we show analytically that that their canonical energy can be made 
negative; the fact that they are unstable (or marginally stable) follows analytically from 
this and numerically from Andersson’s frequency. 

Formal work on the stability of rotating stars does not distinguish between axial and 
polar modes, but explicit studies of nonaxisymmetric instability of stellar modes had been 
restricted to polar perturbations, for which the onset of instablity is very different. There 
are unstable polar modes with arbitrarily small rotation, but as —> 0, the minimum m for 
which a mode is unstable grows without bound. 

In Sect. |2;. we consider the stability of axial perturbations in the Newtonian limit, 
showing that for any nonzero the canonical energy can be made negative for each 
m > 2. The perturbations we consider have angular behavior e im '^ and parity (—l) m+1 , 
corresponding to an l = m mode of a corresponding sperical star. In Sect. [|, we extend the 
computation to general relativity, again finding initial data with negative canonical energy. 
Finally, in Sect. ff|. we estimate the timescale and briefly discuss implications and future 
problems. 

Spacetime indices will be lower-case Greek letters, spatial indices lower-case Latin. We 
will work in gravitational units, setting G — c — 1. 


2. Newtonian stars 


We consider a uniformly rotating, axisymmetric, self-gravitating perfect fluid, described 
by a gravitational potential z/, density p, pressure p and 3-velocity 


v a = 


( 1 ) 


These satisfy the equilibrium equations 


V 2 z/ = 47rp 


( 2 ) 


Va{h - ^v 2 + u) = 0, 


( 3 ) 


together with an equation of state of the form 

p = p{p). 

Here h is the specific enthalpy in a comoving frame, 


( 4 ) 



( 5 ) 
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We will work in the Lagrangian formalism (Friedman and Schutz 1978a) describing 
perturbations in terms of a Lagrangian displacement £ a that connects fluid elements in the 
equilibrium and perturbed star. The Eulerian change SQ in a quantity Q is related to its 
Lagrangian change A Q by 

A Q = 5Q + (6) 

where is the Lie derivative along £ a . 

The displacement £ a determines the fluid perturbation: 


Av a = d t C, 

A = A = _ v „ r , 

7 P P 

For angular dependence we have the corresponding Eulerian changes 


5v a 

5p 

dp 


(d t + imtyC 

-Va(pC) 
dp 


dp 


dp ; 


and the change in the potential is given by 


V 2 5u = 47 t Sp. 


(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 


( 12 ) 


An axial perturbation of a spherical Newtonian star is a change in the fluid’s velocity 
of the form 

Sv a = C (r) e abc V b rW c Yr. (13) 

That is, it is a vector with arbitrary r-dependence, whose behavior on a symmetry 
two-sphere is given by c ab V b YL 0 

The time independence of the perturbation can be seen as follows. An axial 
displacement vector has the form 


C = ((r,t)e abc V b rVcY l m . 


(14) 


2 One can regard a tensor on the two-sphere as a tensor on R 3 that is is orthogonal in all of its indices to 

V a r and is Lie derived by V a r. The antisymmetric tensor €ab associated with the unit metric cab on the 
two-sphere is then identified with e ab = e 0 {, c V a r, and e AB V bY™ is identified with e ab V{,y'; m . 
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The perturbed scalars (dp, dp, dv) must vanish, because they are proportional to Y™ in an 
(l, m ) representation and thus cannot be axial. Explicitly, from the relations 


v„r = o, ev aP = o, 


(15) 


it follows that the Eulerian change in p —and hence in p and v —vanishes. Then only the 
velocity changes, with 

dv a = d t C- (16) 

and the perturbed Euler equation, 

1 


S[(d, + £„)»„ + V„(fc - 5ji 2 + v)\ = 0, 


(17) 


becomes simply 

d t dv a = 0. (18) 

Thus for a spherical Newtonian star, axial perturbations are time-independent changes of 
the velocity, having the form ([13]). 

We will shortly derive the axial instability in a way that does not rely on normal 
modes. Let us first however note that the known frequency of axial modes suffices to imply 
their instability for slowly rotating models. 

For a slowly rotating star, the frequency cr of these modes is no longer zero. Instead, 
conservation of circulation implies in the Newtonian limit (Papaloizou and Pringle 1978) 
the relation 


a + rriQ = 


2mCl 


(19) 


1(1 + 1 )’ 

where the perturbation’s time dependence is e iat . 

Eq. (0) is a consequence of conservation of circulation, the curl of Eq. & 

(d t + imQ)e abc V b Av c = 0. (20) 

(To obtain this equation, note that in Eq. (|T7|) d can be replaced by A, because the 
equation is satisfied by the equilibrium configuration.) For £ a of the form (p~4|), we obtain 


abc 


V b Av c V a r = (d t + Qd rt> )e rH (d e ^ - d^e) - 211 sin 6g 


= i[(a + mQ)l(l + 1) - 2mfl]4^ m 


Eq. (2T[) then has as a solution 


a = —mVt 


1(1 + 1 ) 


( 21 ) 


( 22 ) 
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The frequency cr is negative when l > 2 while the frequency cr + mil relative to a comoving 
observer is positive (Eq. 0). Because modes that travel backwards relative to the star, are 
pulled forwards relative to an inertial frame, they are formally unstable for arbitrarily small 
values of the star’s rotation. 

Stability is governed by the sign of the canonical energy E c , when expressed in terms 
of canonical displacements £ a (Friedman and Schutz 1978). That is, a model is unstable to 
perturbations with angular dependence e im ^ when there is initial data with this angular 
behavior for which E c (£) < 0. The canonical energy has the form 

E ,„ = i/H«| 2 -Pl<'-V?| 2 + 7J>|V-f| 2 

+ (e • Vpv • i + i • Vpv • i 

+ re b (V a V b p + pV a V 6 u) - ^\V6u\ 2 }dV. (23) 

47T 

We evaluate it to order fi 2 . 

For a slowly rotating star, 

P — Po + 0(f2 2 ), p = po + 0(Q 2 ), v = v o + 0(fl 2 ), (24) 


with po, po and vq values for the corresponding spherical star with the same mass. Then for 
£ a of the form (]T4|), we have 


e-Vp = 0(fi 2 ), e-Vp = 0(fi 2 ), e-Vz/ = 0(fi 2 ), 


(25) 


and 

Ec = -jdV [pier - p\v ■ vei 2 + re 6 (V a V b p + pV Q V Q z/)] + 0(fi 4 ). (26) 

In the final term, 

V 0 V b p + pV a V 6 z/ = pV a (-V 6 p + V b u ) + -V a pV b p 

\P J P 

= ^pV a V 6 U 2 + -V a pV b p, (27) 

z p 

where the equation of hydrostatic equilibrium, (Q), is used to obtain the last equality. Then 

e a £ b (V„V b p + pVaVfti/) = ^ P n 2 cev a v b (zu 2 ), (28) 


with til = r sin 0 the distance from the symmetry axis. 
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After a straightforward computation, we find 


E r =l 


A T/ n 


i A\2 


(29) 


An initial data set (£, £) is canonical if it satisfies 

q a = e abc VbA ^ c = Q ( 30 ) 

Writing v a = Q(p a for the equilibrium star and z a for the unit vector V a z, we obtain 

q a = (d t + fl£ 0 )e a6c V b £ c - (31) 

In a cylindrical chart w, z, (p, the components of q a are 

q l = (d t + - d k ij) - 2fld z C, i,j, k cyclic. (32) 


To show the existence of canonical initial data that makes E c negative for each m 
we specialize to the case l = m. A displacement £ a given by eq. (O) with l = m will be 
canonical if for suitable cr, 

C (r) OC r m+1 , d t ( = icr(. (33) 

We can choose 


r 


iw m ~ l ze irn(l> 


-i w m e im(t> 

-w m - 2 ze im 


(34) 

(35) 

(36) 


We then have 


q™ = i [(a + + 1) — 2mVl]w m 1 e im ^ > 

q z = 0 

q^ = —[(a + mQ)(m + l) — 2mQ\zu m ~ 2 e irn ‘t > , 


(37) 

(38) 

(39) 


and q a = 0, when 


cr = — mO 



2 

m(m + 1) 


(40) 


The agreement with Eq. (|22|) reflects the fact that the displacement vector of a normal 
mode is canonical. 








Armed with £“, we can verify that the canonical energy is negative. From 


have, 


and 


|£| 2 = ^ sm 2 *"*- 1 ) 0 (l + cos 2 0), 
d 2 

r 2 sin 9 cos 9 |«f^| = sin 2 ^" 1 ) 9 cos 2 9. 


we 

(41) 

(42) 


Then, explicitly performing the 9 integration, we find that the canonical energy of Eq. 
is given by 


Ec=~^j dVp\i\‘ 


implying 


(mQ — a ){mVt + a) — 


E c < 0, 


2 m 


2 rn + 1 


-n 2 


(43) 

(44) 


for a satisfying Eq. (|T9|). 


For a normal mode, the canonical condition is equivalent to the curl of the perturbed 
Euler equations. In fact, the displacement £ a given by Eqs. ( fil - flbl ) (or, equivalently, by 
Eqs. (0 and (p3|)) satisfies the full Euler equations and is the correct form of an l = m 
axial mode up to terms of order fl 3 . 


3. Relativistic stars 

The spacetime of a uniformly rotating relativistic star is described, to lowest order in 
angular velocity, by a metric 

ds 2 = —e 2v dt 2 + e 2 ^(d0 — udt) 2 + e 2X dr 2 + e 2 ^d9 2 , (45) 

with commuting Killing vectors t a and (j) a . The fluid has four-velocity 

u a = y t (t a + flcj) a ), (46) 

and its energy density e, pressure p, and specific entropy s, satisfy an equation of state, 

e = e(p, s) (47) 

(neutron stars are effectively isentropic, with no dependence on s). The field equations 

G a /3 = 8irT a g (48) 

imply the equation of hydrostatic equilibrium, 


V a (h + log-?/) = 0, 


(49) 






where 


dp 

e+p 


(50) 


h = 


As in the Newtonian limit, an axial perturbation of a spherical relativistic star involves 
a time-independent change in the fluid’s 4-velocity, and no change in the star’s pressure or 
density. (The metric, however, does change, as will be discussed below in the treatment of 
initial data describing a perturbed rotating star.) The change in the 4-velocity has the form 

= C(r) (*‘V s lV,rV,l': i (51) 

(52) 

(53) 


where 

with 


e a/3 := 


ra = e V a r. 


For a rotating star, the frequency cr of these modes is no longer zero. To show that 
axial perturbations with arbitrary values of m are unstable for slowly rotating perfect-fluid 
models, we shall again consider perturbations corresponding to a Lagrangian displacement 
with initial data of the form 


e = cm 

£,{" = r)f« 


(54) 

(55) 


The dependence of cr on r arises because it is easier to work with initial data that is not 
data for a pure outgoing mode. 


Initial data for the perturbed star and geometry is a set (£ a , £ t £ a , h a p, £th a p), where 
h a p = 5g a p, satisfying the initial value equations, 


5(G a0 - 8nT ats )\7 l3 t = 0, (56) 

and the equations governing a canonical displacement. Given such a set, we have (see 
Friedman and Ipser 1992 for a review of the perturbation formalism used here) 


^9a/3 hap + 2 V 3), 


and 


Au a = -Ab/V'A. 9;37 , 

— = Al = — = -W 

7 p e+p n 2 


(57) 

(58) 


n 


(59) 
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A spherical star has metric 

ds 2 = —e 2l/ dt 2 + e 2X dr 2 + r 2 (dd 2 + sin 2 6d(j) 2 ) (60) 

and four-velocity 

u a = e~ u t a . (61) 

Axial perturbations of the metric can have nonzero vector contributions 

h a pt l3 i h aj3 r 0 (62) 

proportional to the vector 

(63) 

and a nonzero tensor contribution proportional to 

^v / 3 ) v 7 rr- (64) 

The form of implies 

A u a = 0, (65) 

5u a = -e~ u £ t C, ( 66 ) 

and 

5e = 5p = 0. (67) 


A slowly rotating star has metric for which u is O(Q) and the remaining departure 
from spherical is 0 (fl 2 ): 

ds 2 = —e 2v {l + 2 h)dt 2 + r 2 (l + 2k) sin 2 6(d(j) — udt) 2 + e 2A (l + 2 £)dr 2 + r 2 ( 1 + 2 k)d6 2 ] (68) 


we will need only the fact that h, k,£ are 0(fl 2 ), not their explicit values (Hartle 1967). For 
a displacement of the form (|54p, we have 


a u a = o(n 2 ) 

Se = o(Q 2 ), Sp = o(n 2 ) 


(69) 

(70) 


The condition governing a canonical displacement is again related to conservation of 
circulation, which has the relativistic form (Friedman 1978) 


£\i^a/3 0 ) 


( 71 ) 
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where uj a p is the relativistic vorticity, 

^a0 2V^ 


e + p 


n 


-up] 


Then canonical initial data satisfies 


Now 


with 


AcUq;^ 2V[ a 


0 

e+p 


n 


a Uf3] ) + o(n 2 ), 


A u a = A (g a pu 0 ) = Agapu 13 + 0(fi 2 ) 

= 2V {a ^ + 0(n 2 ). 


t 41 = - sin m_2 9 cos 9 

£t£ = ia£. 


From Eq. ([f4D, we obtain 


satished when 


a + mfl = 


2 m (Q — a;) 


(72) 

(73) 

(74) 


(75) 


As in the Newtonian case, we will use data with l = m, writing £ Q in the form 

5" = i sin ™- 1 0 e im * ++ (76) 


(77) 

(78) 


Au eij> = 25 [( j +^-A u+ = 0 : 

u 1 sin m+1 9 e 1 ™^ + l)(cr + m£l) + 2 im(Q — a;)] = 0, (79) 


(80) 


m(m + 1) 

Thus, as noted earlier, we have not supplied data for a normal mode, but data for which cr 
is a function of r, with 


d r cr = 


-d r oj. 


m + 1 

The remaining components of the canonical condition ( |73| ) are 

Auj r< j, = 0 : 


(81) 


d r log ( " ] r 'j AUrj, + d r ( AUrf) - (Au r ) = 0, 


(82) 
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and 


/\cu r g — 0 : 

d r log Aw*? + (Atte) — dg (A u r ) = 0. (83) 

Using the equation of hydrostatic equilibrium in the form 


one obtains from the remaining components of the canonical condition a single equation for 

c, 


drC -—;—K 

m + 1 

with solution, 


m + 1 


d r u - -d r {r 2 (VL — u>)) — 2d r u 


'2 (n-u) 

m + 1 


— (Q — u) 


c = (ft - cn)— 1+(m+1) / 2 r m+ i exp [-(m - l)i/]. 


= 0, (84) 


(85) 


We have obtained a canonical displacement, but must still satisfy the initial value 
equations. We will do so by choosing a metric perturbation for which h a p = 0, £ t h a p ^ 0. 
The initial value equations can be obtained from the action for the perturbation equations, 
which we will use in any event to write down the canonical energy: 

I = J d 4 xy/=g£(£, h f, h) ( 86 ) 

where 

C{lh\Z,h) = U^V q ^V 7 6 + ^ 75 (VV 7 6 + ^V 7 6) 

527T 

-T^R^sCZ 5 + Qw a ^ s - -^-G Q ^ S ) h a phys 
-^V 7 T^(Vf + Vn. (87) 

with 

G 1 *' = + 1(2+ 2IV s g a ® - 3R a( ' r g‘ >13 - 3 R^g**) 

+ + ( 88 ) 

f/o/Sli _ ( e + p )u<*u>ql ,s + p(jV - g“ s g l, ~') - 7P rty 


(89) 
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and 


2V a W = (e + p){u a vkq p& + vPvkq 


ol 6 


u a u t> q< s ) - -IP rtV 


(90) 


The perturbed initial value equations are given by 


0 


51 


5k 


Vflt. 


a/3 


(91) 


When h a p vanishes on the intial data surface, the Hamiltonian constraint vanishes. We 
choose h a g to have no vector contribution, writing 

Kp = h(r )e (o 7 V p)\' 1 Y l m . (92) 


The momentum constraint then has, to 0(H), the form 


(e + pK£ u r - e ta ^\Wph lS = 0, 

OZ7T 


or 


2 i/ 


Cl e 

i(a + mO)(e + p)(u t ) 2 — + —— (m + l)(m + 2 )——h = 0. 

r z ill r r 4 


Thus the initial value equations are satisfied by choosing 


h(r) = —32m 


a + mfl 2 


(l - 1 )(l + 2 ) 


r (e+p)(. 


(93) 

(94) 


(95) 


We can now compute the canonical energy, whose relativistic form is 

E c = f dS e [U&%V^s + V**h 1 sZ p -^e^XeW\hi S Vph m --f£]. 
Js Hit 


For a displacement having the axial form (154|) , we obtain the simpler expression 
E c = JdVe v [U t ^ s £ t ^W^ s -^e t ^e t6 \£ t h lS £ t h eri 

- + \T aP R aiP5 Ce] 

1 + £f2 + Tfz + £g ], 


(96) 


(97) 


where the terms Efi indicate contributions from the fluid, £ g the contribution from the 
perturbed metric. (This split into fluid and metric is gauge-dependent). Evaluating each 
term, we obtain 


£fi = C/ w £.&V 7 6 = -ia(e + p)u%u<V i e + 0(S! 4 ), 


(98) 










(99) 


Sf 2 = --[/^V q ^V 7 6 = 


^e + p)\u-VZ \ 2 + ^pVj g V^ a 

+ 0{n 4 ). 


£ f3 = \T^R aiP5 Ce = - i(e + p)uN^u 5 V^s - 

+ l -{e + p) (/V a ^ 7 V 7 « Q + m q V q ^V 7 /) 

+ e- v D a (e v T>H®Vfa). 


( 100 ) 


In this final equality, D a is the derivative operator with respect to the 3-metric on the 
t ^constant initial value surface. Then 

£f = £fi J r£f2 J r£ 


= ^(e + p^V (a £ g) £ u C - ia(e + p)u^ /3 u 7 V 7 ^ /3 


+e~ l 'D a A a , 


^(cr - m£l)(e + p^t^pvAV ^ + ^(cr + m£l)(e + p)u^“^ /3 V a u g 


where 


1 V^ 7 + 2 (o' + mfi)(e + ?)«*£“«%]. 


The remaining term is 

*. = 1 


J^e te \ EthyMhsr, = 7[7~ e ~ 2v ^ e 75 £ t fi a7 £ t h gs , 
bin 64n 


where e a p is the metric on a 2-sphere of radius r, 

e a/3 = r 2 (VJV g 6 + sin 2 d V a (/> V/ 30 ). 

Thus 


E c = 

where the terms have opposite signs: 


+ £?); 


£f = — — {mQ, — cr)(mQ + er)e ^(e + p)|£| 2 

+2mQ(n — cu)e _2!y (e + p)r 2 sin 0| cos0£ 0 £^| < 0 

4 = ^e- 2v e^e / % s h cer > 0. 


( 101 ) 

( 102 ) 

(103) 

(104) 

(105) 


(106) 

(107) 
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The fluid contribution Sf has the same form as the canonical energy 
modes, from which the sign of (|106|) follows. 


for Newtonian 


The sign of E c is negative because the fluid contribution is dominant. Writing 


S f = 


m + 1 


(fi — u) 


mil 


— u 
m + 1 


cos 2 9 ( m 2 Vt + 


— u 
m + 1). 


x 


c 2 

3 2u (e + p)-^ sin 




S g = 327r 


m 


2m—2 n 
1 2 


(108) 


m + l)(m + 2) 
e~ 2u (e + p) 2 ( 2 sin 2m 9 1 1 + 


( fi - u ;) 2 
8 cos 2 9 ' 


we have 


Idns g 


16 


sin 4 0 


27r(e + p)r 2 < 27r(e + p)r 2 . 


(109) 


( 110 ) 


\JdQSf\ \{m — l)(m + 2)] 2 

The maximum value of 27r(e + p)r 2 , for any equation of state, occurs for the maximum mass 
star. We directly calculate 27r(e + p)r 2 for the maximum mass star using equations of state 
spanning a wide range of stiffness. For equations of state G (excessively soft), C and L 
(excessively stiff) from the Arnett and Bowers (1977) catalogue we find 


27r(e + p)r 2 < 0.4, 


for all radius, implying 


E c < 0. 


(111) 

( 112 ) 


4. Discussion 

Although a mode with negative canonical energy is unstable, the growth times of axial 
modes appear to be proportional to a high power of for slowly rotating stars. Because 
the nonaxisymmetric instability of rotating perfect-fluid models is driven by gravitational 
radiation, its growth time for modes with a given m is ordinarily shortest for modes that 
correspond to the smallest value of / for a given m, namely l = \m\. For axial perturbations 
of slowly rotating stars, however, modes corresponding to a given l and m for the spherical 
star give rise to density perturbations associated with Y^. For nonisentropic stars, axial 
modes with l = \m\ + 1 give rise to Y™ density perturbations, and as a result, these 
modes have the shortest growth time. However, if the star is isentropic, it can be shown 
(Provost et. al. (1981)) that the only axial normal modes of a Newtonian star have 
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l = \m\, which contribute a m + 1 multipole to the density perturbation. As a result, 
the gravitational radiation from mass multipoles is dominated by radiation from current 
multipoles (Lindblom, Owen & Morsink 1998). 

We provide here a rough estimate of the time scale for the axial instability of a mode 
with angular dependence exp (im<j>), for a slowly rotating, nearly Newtonian model. A mode 
characterized by a displacement vector £ has energy of order 

E -M0 2 £ 2 . (113) 


The mode is driven by gravitational radiation, with power dominated by the lowest order 
current multipole of order 

f ~-(!>) 2 aw) 

where J” 1 is the lowest current multipole perturbation of the star. We have 

~ J p Sv a e a bc V b rV c Y l m r l dV. (115) 


With our choice of £, 5v a is 0(h2)Q 

jr ~ f ^R l+ 2 


(116) 


Finally, the growth time is of order 


dE Mfl 2 £ 2 
' dt ~ n 2m +' 2 \Jr\ 2 
M 

(m) 2m + 2 ' 


(117) 


The precise calculation is given by Lindblom et. al. (1998). The timescale appears to have 
the same hi dependence for relativistic stars (Friedman, Lockitch & Morsink 1998). 

The imaginary part of the frequency that Andersson (1997) computes is 0(h2 2 ), and 
our timescale estimate suggests that its magnitude is an artifact of an 0(fl 2 ) approximation 
method. But the sign of the imaginary part that his method obtains may correctly diagnose 
instability of a mode in the slow-rotation approximation. 


3 It would be more natural from a physical standpoint to take Sv a = 0(1), £ = 0(Q 3 ) because as 17 —> 0, 
the perturbation goes over to a nonvanishing, time-independent perturbation of the spherical star. 
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The discussion so far has ignored viscosity, and, as in the case of polar modes, a long 
growth time for slow rotation implies that viscosity will enforce stability except for hot, 
rapidly rotating neutron stars (see Lindblom and Mendell 1995 and references therein). In 
particular, because the time-independent axial modes of a spherical star have nonvanishing 
shear, viscous damping must damp the instability for slow rotation. 

For rapidly rotating stars, however, instability points have been computed numerically 
only for polar modes, initially for Newtonian models (Managan 1985, Imamura et al. 
1985) and more recently for relativistic models in a Cowling approximation by Yoshida 
and Eriguchi (1997) and in the full theory by Stergioulas and Friedman (Stergioulas 1996, 
Stergioulas and Friedman 1997). Estimates of instability points when realistic values 
viscosity is included have similarly been restricted to polar perturbations. Extending these 
estimates to axial perturbations should then be straightforward. But for both axial and 
polar modes, a precise computation of instability points when viscosity is included is a more 
difficult task, apparently requiring the construction of outgoing modes and determination 
of their complex frequencies. 


We are grateful to Nils Andersson for acquainting us with his results and for several 
helpful discussions. This work was supported in part by NSF Grant No. PHY95-07740 and 
by NSERC of Canada. 
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